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ABSTRACT 


In this thesis a new antenna configuration; an array of 
electrically small dipoles, is investigated as a possible solution 

J. 

to the requirement o£ a broad band antenna with a dipole like 
pattern. Multi octave bandwidth is what the proposed antenna 
configuration aims at. The configuration consists of a linear- 
array of electrically small dipoles fed by a specially designed 
feed network, to maintain an input match as well as a dipole like 
current distribution over the entire band of frequencies . 

The analysis is carried out in three parts. First, an 
individual electrically small dipole is analyzed to obtain Its 
input characteristics as a function of frequency. In the second 
part mutual coupling between such dipoles kept close edge wise is 
analyzed. This resulted in characterizing the entire antenna array 
in terms of n-port scattering matrix, with each element of the 
matrix being a function of frequency. The third part consists of 
designing the terminal parameters of the feed network which 
provides the necessary driving point Impedances and excitation 
voltages to the antenna element, which in turn maintains the 
desired current distribution on the antenna array. This analysis 
involves solution of simultaneous coupled equations which are 
solved using an iterative procedure. The best possible match at 
the feed port of the feed network as well as best match to the 



desired driving point impedances at the antenna terminals, are 
obtained using an optimization algorithm based on a pattern search 
t echnique . 

The analysis given in this thesis provides complete 
procedure for calculating optimum terminal parameters or in other 
words, the complete scatter matrix of the feed network as well as 
the antenna array. The actual design of the feed network Is beyond 
the scope of this thesis. 



CHAPTER -1 INTRODUCTION. 


1.1 Introduction and literature survey : 

Direction finding systems C DFD require omni di recti onal pattern 
in azimuthal direction and preferably constant beam width in 
elevation over a wide band of frequency. These systems normally 
use resistive wire dipoles or bi coni cal antennas for wide band 
coverage because a normal resonant X-'2 dipole has only 15% to 20% 
bandwidth. One of the drawbacks of resistive dipoles is that they 
are inefficient. The work reported in this thesis is an effort 
towards obtaining a highly efficient antenna while retaining the 
properties of wide band * omni di recti onal coverage pattern etc. * 
required for a DF system. In addition to these requirements 
desirable property of retaining the beam width constant in 
elevation over the band of frequencies of operation* has also 
been taken into consi der ati on in this study. 

Literature is abundant in the area of dipole analysis. 
Effort towards broad banding a dipole has been reported by several 
authors. Only the ones relevant to our study have been reviewed 
here. King and Wu Cl] have shown that* dipoles will have wider 
bandwidth if the standing wave type of current distribution is 
converted into travelling wave type. This is due to the fact that* 
input impedance is a strong function of frequency for standing 



wave current, distribution* whereas travelling wave current makes 
it less dependant on frequency. 

L. C. £hen [S3 suggested a resistive loading technique to 
extend the usable bandwidth of dipoles. Obvious disadvantage with 
resistive loading technique is reduced efficiency. A study of 
reactive loading technique was carried out by B. L. Rao [33. Most- 
useful antenna * reported by Rao is an Exponentially Tapered 
Reaetively Loaded Antenna [ EXTRA] . This is reported to have 3:1 
band wi d t h . 

B. J. Strait and K.Hiraawsa [43 have investigated multiple fed 
long wire antennas. They could achieve arbitrary high directivity 
in some direction by properly choosing the current distribution. 
Tliis type of arbitrary current distribution needs optimized 
excitations and source impedances. Recently* SaudyL B3 has 

employed the came technique in O.SX* IX* 1 . SX and 3X dipoles to 
optimize gain and bandwidth. His analysis shows that an 

improvement is possible only at the expense of a complicated feed 
network and increased side lobe level. 



1.2 Problem definition and proposed approach to the solution : 

The aim of this study is to develop a radiator with dipole 
like pattern and having a multi octave bandwidth. The requirement 
in DF systems is that the coverage be omnidirectional in azimuth 
and preferably constant beamwidth in elevation. A typical \/2 
dipole pattern fits the pattern criteria very well except that 
dipole bandwidth is very small. 

It is a well known that for good radiating efficiency the 
size of the antenna has to be at least comparable to quarter 
wavelength. A half wave dipole gives a good efficiency with about 
15% to 20% bandwidth. Because of the resonant nature of the dipole 
fed at one point the bandwidth is limited. There have been studies 
in which a dipole is made to radiate efficiently at two or more 
frequency bands by inserting filter networks in the dipoles such 
that the dipole is resonant at different frequencies. Essentially 
these filter networks cut off the end sections at a higher 
frequency confining the current distribution to only the center 
sections. The work reported in this thesis is an extension of this 
idea in some way. 

Our requirement here is to get a continuous broadband 
operation while maintaining dipole like pattern. Two conditions to 
be satisfied by the antenna become clear from this. One is that 



the current distribution should be similar to a dipole current 
distribution at all frequencies in the operating band > which 
essentially means that we have to restrict the effective radiating 
length of the antenna to at all the frequency. Thi s also 

ensures high radiation efficiency for the antenna. The second 
condition is that the input should be matched to the source over 
the entire operating frequency band. Thi s is sought to be 

achieved by the proposed antenna configuration. 

The proposed configuration is that a linear antenna is built 
out of small dipoles* each with its feed point at the center of 
the dipole. The total length of the array of these dipoles is 
made equal to half wavelength at the lowest frequency of 
operation* so that it remains an efficient radiator over the 

entire band. Obviously* the individual dipole elements are 

electrically small at the lowest frequency* and would have very 
low efficiency. However* when all the elements taken together 
would produce a current distribution which resembles X^“S resonant 
dipole current distribution* and would be efficient as an array. 
To obtain a match to the source over the frequency band of 

interest* a 1 eed network with specially designed characteri sti cs 
is inserted between antenna elements and the source. This 

essentially would be a l:n power divider network with the power- 
distribution ratio of each port having a specially desiqned 



function of frequency which in turn maintains the antenna current 
distribution to conform to a dipole like distribution with 
frequency. As frequency is increased from the lower end of the 
operating band* the extreme elements are progressively discon 
nee ted or the currents on these extreme elements are forced to 
zero* to keep the effective radiating length equal to X.^-2. Having 
defined the configuration capable of giving the desired 

e har- ae t er- istics of dipole like pattern and broad bandwidth* the 
next- step is to evolve a method to design such an antenna. 

The configuration we have is a linear array of electrically 

small dipoles with close spacings- and is fed by a feed network 

with some special characteristics. The problem is one of 

determining the length of the array* size of the elemental 
dipoles* number of dipoles* their spacings and the feed network 

parameters* given an operating frequency band. First two 

parameters can be easily arrived at without much anal ysi s. from the 
current distribution requirement the total length of the array is 
fixed to be at lowest operating frequency. At the hi ghest 

oper-ati ng frequency* it is envisaged that only one element will 
radiate* thus* the center element length is fixed to be at the 

highest frequency. The same length is also taken to be element 
length. It is not necessary that this be so one can have elements 
of different sizes making up the array. It is not yet determined 



what, is -the? best, or optimum way of choosing the element lengths. 
In this study we have assumed identical elements for simplicity of 
analysis. Further spacing is assumed to be slightly longer than 
the element lengths so that t these elements are end to end* but 
do not touch each other. This gives nearly continuous current 
distribution over the array length. This* of course* imposes a 
condition that the array length be equal to integral multiples of 
at the highest frequency. One can chose the highest and 

lowest frequencies such that this criteria is satisfied. Once 
element length and array length is determined* the number of 

elements is also automatically fixed. The number of elements is 
made odd so that we have a center element which alone radiates at 
the highest frequency of operation. 

The analysis required to arrive at the design of the entire 
antenna system can be divided into several parts. Because the 

elements are electrically small and are closely spaced » the effect 
of mutual coupling is very strong* and it is important to 
characterize this accurately. The analysis is done in the 

f ol 1 owi ng stages . 

al> Analysis of electrically small individual* isolated radiator 

for its input impedance behavior over the frequency band of 

interest. 

bl) Analysis of mutual coupling between two such elements. 



cD Input impedance behavior with frequency of "the elements in an 
array. 

dZ> Calculation of driving point impedances and excitation voltages 
required at each of the element inputs to sustain the desired 
dipole like current distribution. 

el) Deter mi nation of feed network parameters which will provide the 
necessary excitation voltages and driving point impedances as a 
function of frequency 

f2> Final ly* the design of the feed network itself from the 
required terminal parameters as arrived at step<e3 . 

These steps Ca2> through C&J have been formulated and 
discussed in this thesis. The last step CfO has not been 
attempted. 


1.3 Present investigations s 

This work on array of electrically small radiators considers 
dipoles as individual radiators. To start with* electrically 
small dipole is analysed in chapter -2. The input impedance is 
evaluated as a function of frequency* using well known Induced EMF 
method. Results of this chapter show that* electrically small 
radiators are highly reaotive and are not efficient if operated as 



an individual element antenna. The efficiency increases as the 
length of the element becomes comparable to resonant length X.y'S. 

Chapter -3 describes about mutual coupling aspects of closely 
spaced electrically small dipoles. Here also Induced EMF method 
is used to calculate mutual impedances. Then the equivalent 
circuit for the array is developed to evaluate driving point 
impedances and to study the effect of some termination on coupled 
antennas . 

Mext chapter is entirely devoted to feeding aspects of 
closely spaced array. Feed network parameters are determined such 
that a single source excites all the elements of array. Feed 

network synthesis uses an optimisation routine to achieve minimum 
VSWR at the primary feed point. 

This chapter is followed by a chapter on summary and 
conclusions. The synthesized S-matrix in conjugation with antenna 
elements are evaluated at many spot frequencies throughout the 
frequency of interest. Radiation pattern at different frequencies 
are also presented for the designed antenna. Some suggestions for 
further work are included at the end of the chapter. 



CHAPTER 2.INPUT IMPEDANCE OF SMALL DlPOLES 


2. 1 Introduciion : 

An antenna can be represented by an equivalent Impedance at 
the input terminals of the antenna, as far as the feed network Is 
concerned . Impedance of an antenna at a single frequency can be 

written as: 

Z = ( R + R ) + j X ... (2.1) 

•cl r L a 

where , 

R - Radiation resistance of the antenna 
R^ - Loss resistance of the antenna 
X - Reactance of the antenna 

CL 

This impedance 'Z ’at its terminals is called "Self 

CL 

impedance”, which is the input impedance of the antenna in an 
isolated free space. In the presence of any interfering element 
(say some other antenna or electromagnet i cal ly sensitive material) 
the input impedance is modified by them and the impedance 
seen at its terminals is known to be "Driving Point Impedance”. 
There are various methods to calculate the self Impedance of the 
antenna. Generally these methods can be catogorised Into 



1. Boundary value method. 

2. Equivalent Transmission Line method. 

3. Poynting Vector Method. 

4. Numerical Methods. 

Out of these Poynting Vector Method [6] is the most popular- 
one because of its simplicity. In the Poynting Vector Method, 
Power density is integrated, over a closed surface enclosing 
the antenna to yield the total power. Power divided by the square 
of the assumed current at the input terminals gives the 
input impedance. 

2.2 Near- f ields of' a dipole s 

Referring to the figure 2.1 and assuming a sinusoidal current 
distribution along the element, as given by 

A 

I (x'.y' ,zf) = a I sin [ k (1/2 - * )] 0 < z'< ■ 1/2 

91 O 

Aw 

= a Io sin [ k (1/2 + z' )] -1/2 < z’< 0 ; 

z 

. . .( 2 . 2 ) 

Vector potential ’A’ can be written as, 

■+ A 

A = a A 

z z 



X 
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PCOCyV, 2 *) 



Fig 2.1 Geometry of DIPOLE fox- near field calculations. 
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If (X',Y’,Z’) are the source co-ordinates. 


r = J (x 2 + y 2 + (z-z' ) 2 ) 

= VI? + (z-z’) 2 

Magnetic field H is related to vector potential *A* 

H = 1/aj V x A 


For an A component with no $ variations, in 

sc 

co-ordinates, H becomes. 



dz’ 


. . .( 2 . 3 ) 


. . . ( 2 . 4 ) 


by 


. . . (2.5) 

cylindrical 


. . . ( 2 . 6 ) 



Since the field is not a function of 4 > , because of azimuthal 


symmetry, the observations can be made at any angle 4 > . 
phi is equal to n / 2, x,y becomes, 

X = p Cos = 0 ; 

Y = p Sin - p ; and S/Sp = S/Sy. 
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After some mathematical manipulation, H is simplified as 


H 
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2 cos (kl/2) e 


- . -( 


Where R ,R ,R are shown in the fig. 2 . 1 . Corresponding 
1 2 

field components are given by, £ = V x n / 


* 

E 

P 
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a h , 
4 > 

* y 


Choosing 


dz’ 


. (2.7) 



. 8 ) 

electric 

C2.9) 

.( 2 . 10 ) 


J £ 




These are the required expressions for the fields radiated by a 
dipole . 

2.3 Input Impedance by Induced EMF method : 

Induced EMF method [7] can be briefly explained as follows. An 
EMF is applied at the terminals of the antenna, which produces a 



current . The current produces an electric field E which 

in turn induces a field E - at the conductor surface such that 

the boundary conditions are satisfied. For a perfect conductor 

E.=E+E,=0 at the surface 

zt z zi 

Therefore E = -E , 
z zi 

Potential difference dV^ over an incremental length 

dV = -E dz’ 
z z 

to the current maximum by the 


. . .( 2 . 16 ) 


This voltage is related 
transfer impedance i.e., 

dV 

-7 * 


. . .( 2 . 14 ) 

dz ' is 
...( 2 . 15 ) 


Because of reciprocity, Z = V /I 

*- m z 

voltage at current maximum. 

I dV = V dl 

* a mm 


where, V is the 

m 


. . .( 2 . 17 ) 


Radiation impedance referred to current maximum is defined as, 


V 

dV 

m 

m 

I 

dl 

m 

m 


...( 2 . 18 ) 



From equation (2.17) and equation (2.18) 


I 


dV = — dV 

m j z 

01 


(2.19) 


L S2. 
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in 
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dV 


m 




1x^2 


- / 


I E dz’ 
z z 


.(2.20) 


in 


-Vya 


Therefore, 


0t 


1 

[* 

m 


l y2 


-LX2 


I E dz’ 
z z 


( 2 . 21 ) 


where , 


I sin [ k (1/2 - tz'l )] 
m 


. . ( 2 . 22 ) 


and E is E field produced by the above current flowing in a 

m 

{ i 1 anient placed along the axis of the wire, (derived in last 
section) . 


Substituting 
manipulating, real 


for E and I in equation (2.21) 

5E 2 

and imaginary parts of impedance can 


and 

be 


expressed as: 



= / C + Ln (kl) - Ci (kl) + 


1/2 Sin (kl)£ Si (2kl) - 2 Si(kl) j + 

1/2 Cos (kl)[ C + In (kl/2 ) +Ci (2kl) - 2 Ci(kl) 


. . .(2.23) 


2 Si (kl) + 


Cos (kl) ^ 2 Si (kl) - Si ( 2kl ) j 
-Sin (kl) [ 2 Ci (kl) - Ci(2kl) - Ci(2ka Z /D 


where C = 0.5772 (Eulers' constant) 

C and S are Cosine and Sine integrals given by: 
i. i. 


cm - - J 


Cos y 


Cos y 


-(2.24) 


. (2.25) 


S. (X) = 
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Sin y 


. . . (2.26) 
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2. 4, Results and Discussions: 


A program ’Self. Pas’ Is written to calculate the input 
impedance of small dipoles, using equations (2.23) and (2.24). 
Fig (2.2) shows the plot of radiation resistance and input 
reactance as a function of length of dipole. It is observed that 
input impedance is dominated by the reactive part (X » R ) and 

<s. r 

its slope is high. For (l/>0 less than 0.25 antenna is an 


inefficient radiator because R « X . 



Chapter - 3 


Mutual Coupling Between Closely Spaced Electrically Small Dipoles. 


3.1 Introduction: 

Antenna arrays are constructed by arranging N-individual 
elements in some specified configuration. In general, the 
behavior of each of these antennas in the presence of others will 
not be same as its isolated behavior. An antenna which is matched 
when isolated may no longer be so, when placed in an array and 
also the radiation pattern of an element in array will not in 
general be same as the radiation pattern of the element in 
isolated configuration. 


These 

two properties 

arise 

due to 

the 

mutual 

coupling 

between the 

elements of the 

array . 

A study 

of 

mutual 

coupling 


between elements is essential to understand the behavior of the 
array. If the elements are placed close to each other, the effect 
of strong mutual coupling plays major role in array design. 



3.2 Current Distribution In Coupled Dipoles: 


Mutual coupling between the elements in a microwave network can 
be expressed in terms of S-paramet ers or 2-parameters. The 


advantage 

with 

the S-parameter is that it is 

easy 

to 

measure in 

practice 

and 

the 2-parameter is mostly 

used 

in 

theoretical 


computations . Mutual Impedance (Z ) is the ratio of voltages at 

'-i 

the i l port due to the current in the j th port and due to 
reciprocity, 2=2. 

vj jv 

The mutual impedance between two dipoles is calculated with 
the assumption that, current distribution in the first dipole is 
same as the current distribution in an isolated dipole. This is 
true only when the scattered fields due to the second antenna is 
zero, which is not possible. Hence to accurately find out the 
mutual coupling it is necessary to solve for the actual currents 
starting from boundary conditions and Maxwell’s equations. 


King et al [8] have given a rigorous formulation to find out 
the current distribution in coupled dipoles. the analysis leads 
to current as a function of position of antennas, given by a three 
term approximation. 


i 

«n 


C2 J 


A 
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ozu 


+ B 


ozu 


D H 


OZM 


. . . (3.1) 



where are complex constants dependant on the geometry. 

. . . ( 3 . 2 . a) 

... (3.2. b) 

• . . (3 . 2 . c) 

Where m is the primary or coupled antenna, Z represents the 

position along the length of the antenna ’m", h is the length of 

m 

the antenna ' m’ . 

Evaluation of these three components individually, for 
various lengths of the antennas show that the contribution due to 
the last two terms F and H are small, in general, for X./2 dipole 
and negligibly small for dipole lengths << X./2. Hence in the 
present analysis, the current distribution in the dipoles are 
assumed to be of first term ( M ) only. 

3.3 Equivalent Circuit of* Coupled Antennas: 

As microwave engineers feel easy to visualize the terminal 
behaviour of the antenna in terms of equivalent circuit, it is 
desirable to represent coupled antennas also, by means of a 
multiport network. For simplicity two coupled antennae can be 
represented by a two port network. 


the functions M,F,H are given by, 

K = Sin ft ( h. - I Z I ) 

F = Cos ft Z - Cos ft h 

OSZM O m f o m 

H = Cos ft Z 12 - Cos ft h / 2 

OZM O Trt o m 



The elements of T-equlvalent 


circuit are represented In 


a5 


Z-parameters . 


The system of equations, characterizing It is, 



Z 

11 
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21 



. . . (3.3) 


The input Impedance or driving point impedance at port-1 is given 
by 


Z 

INI 



z 

11 


+ 



- - - (3.4) 


In a similar way n-element array can be represented by n-port 
equivalent circuit, characterized by means of self and mutual 
Impedances. For an n-element array coupling matrix can be written 
as 


t V] = [ Z] [ I] ...(3.5) 

T 

[ V] = [ Vi ,Vz ] 

T 

[I] = [ II , 12 , 13 b 4 ] 


where 
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Z^ ^ la the coupled or transfer or mutual impedance between 
porta I to j. Due to the reciprocity principle the coupling 
matrix will be aymmetric about the main diagonal. Thia property ia 
made use of in the design of feeding network in chapter 4. 

The Z-parameters Z ,Z can be computed theoretically as 

V V u j 

described in the following section. 


3.4 Mutual coupling in coll near configuration! 

Mutual coupling in a linear planar array configuration has 
been analyzed extensively as early as 1948 [9]. Generalized 
treatment making use of induced EMF method explained in section 
2.3 goes in the following manner. 




The mutual impedance between the dipoles of fig (3.1) is defined by 


Z 

Si 


V 

24 


I 4 (0) 


. . . (3.7) 


where is the open circuit voltage at the terminals of antenna-2 
due to 1^(0), and 1^(0) is the current at the center of antenna-1. 
The electric field due to the current 1^(0) is known (derived in 
section 2.2 -Near fields of dipole). The parallel component 
along the axis of antenna-2 is calculated with the antenna-2 
removed (temporarily). This gives the open circuit voltage across 
the terminals of antenna-2 as 

H I (z) dz ... (3.8) 

£524 2 

<2Lr»t«fr»ri€i 2 
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So that the mutual impedance becomes, 
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I 4 (0)I a (0) J 
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E I (z) dz 

Z24 2 


. . .(3.9) 


I (z) is the current distribution along antenna-2 given by 
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■«<*> - Si " »[-T - M ) 


. . . (3.10) 


substituting this in the expression for E equation (3.7) becomes 
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... (3.11) 

For the col inear configuration, fig (3.2) real and imaginary P 
of mutual impedances are expressed as 
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Fig 3. a. Coll near configuration. 







where V 

0 
V 
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Kh 

2CK (h+1) ) 


V 2 = 2(K (h- 1) ) 


V 

s 


c h 2 -i 2 )/h 
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Cl(x) and Si(x) are Cosine and Sine integrals given by 
equation (2.25) and equation (2.26). 

3.5 Driving point, impedance: 

Equivalent circuit of two coupled antennae is described in 
section 3.3 and also in fig. 3. 3. This helps to visualize the 
effect of source and load impedances on the behaviour of the 
network. The driving point impedance of the coupled antennae 
should be matched to the source impedance at the respective 
terminals. For a two port network with source (or load) impedance 
Z, L and Z^at portl and port2, excitations V t and V z at portl and 
port2 the driving point impedance can be calculated from the 
equivalent circuit as : 


Z = Z + Z (I /I ) 

id ii m 2 i 


(3.14) 



2 



z 

m 

Z + Z 

22 2L 



---(3-15) 


Driving point impedance at port 2 is given by a similar- 
expression (Interchanging subscripts 1 by 2). This says that 
driving point impedance at portl is a function of load (or source 
impedance at the second port, and exciting voltages at portl and 
port2, besides the mutual coupling impedance. Mutual impedance is 
fixed by the geometry of the array, but excitations and load or 
source impedances Z lL , Z aL are independent variables, which can be 
chosen to match the elements. 

For an n-element array, characterized by the system of 


equations (3.5) driving point impedance is given by. 
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. . . (3.1$) 


3.6 Results and Discussion: 

A program "Mutual . Pas” in Turbopascal version 4 implements 

<2, 3 

equations (3.10) and ( 3 . 1)() to calculate the coupling matrix or Z 
matrix of an n-element array o£ dipoles. For the sake of 

simplicity, important results for two coupled antennae are 
presented in fig. 3. 4(a) & 3.4(b). The coupling is significant for 
the closed spaced dipoles and diminishes rapidly with separation. 
That is, in an array of X/2 dipoles only the immediate, adjacent 
antennae will have significant coupling. Uhen the antenna 
dimensions and spacings are small compared to X/2, coupling is 
stronger between several elements situated within a distance of 
X/2 . 


Fig. 3. 5(a), fig. 3. 5(b) show the plot of variation of driving 
point impedance as a function of currents in the antenna. In the 
previous chapter, the individual el ement ( el ectrical ly small 



radiatetl) was found to be highly reactive and was difficult to 
match. A study of fig. 3. 5 reveal the role of the 


currents In 

the antenna( 

and 

implicitly 

the 

loads ) 

to 

vary 

the 

driving point impedance 

suitable 

for 

matching . 

In 

the 

next 

chapter all 

these properties 

are 

used 

to 

synthesize 

the 

feed 


network. 



CHAPTER -4-. Array of electrically small radiators 


4. 1 Introduction: 

Ae mentioned earlier- in chapt er- 1 , the idea o £ selective 
feeding of a linear array of electrically small dipoles placed 
close to each other is utilized to control the current 
distribution on the antenna. Such a linear array sought to be 
designed to operate over a bandwidth ( f to f # ) with efficient 
radiation characteristics and constant beam width. Selectively 
exciting the array, needs a complicated feed network. This 
chapter discusses a method of arriving at the terminal parameters 
of such a feed network. 

Fig 4.1 shows the arrangement of the linear array and the 
feed network. The box named antenna array encloses an array of n 
elements which can be looked at as a n-port network. These 
n-ports are fed by a (n+l)-port feed network. n-ports of feed 
network are connected to n terminals of antenna array and the 
(n+1 ) th port is connected to the source. It Is at this port that 
we have to retain the match over the entire bandwidth, while 
obtaining the desired radiation characteristics. 

A method is suggested here to synthesize the f^eed network 
terminal parameters for a given antenna array which is 
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characterized in terms of a coupling matrix ([S] or [Z]). The 
breakdown of the problem is as follows: 

1. The Z-matrix characterizing the array of closely spaced 
electrically small radiators is determined using Induced EMF 
method as mentioned in section 3.3. 

2. The terminal currents at each frequency are derived from the 
required current distribution, which gives the constant beam width 
pattern. 

3. Given a set of currents at the n terminals of the antenna 
network, the driving point Impedances are determined as explained 
in sect Ion 3.4. 

4. From the terminal currents and the driving point Impedances, 
the excitations necessary to sustain the given current 
distribution are calculated. 

5. The requirements of the feed network is to provide the 
necessary secondary feed excitations (as calculated in step 4) and 
the driving point impedances (calculated In step 3) and to have 
minimum VSUR at the (n+l)th port; the primary feed point. 


This chapter describes how to calculate the reflection 


coefficients at the (n+1 ) th port of feed network, and how to meet 
the requirements mentioned in step 5, while designing the feed 
network . 


4. 2 Ref lection coef f icient in a loaded networks 

Fig 4.2 shows a , a and b ,b as the incident and reflected 

a ft a a 

waves at port 1 and port 2 of a two port network. a’a’and b’.b’ 
are the incident and reflected waves at load/source i . e . , at port 
1 and port 2 . Uhen the loads are connected to the networks we can 
equate 


b 

i 



a 

A 


b * 


Therefore the load reflection coefficient can be written as: 


r 

L, 


b * 
2 


/ a' 


/ b 


. . .(4.1) 


This network can be characterized in terms of a 2x2 S-parameter 
matrix. Uhen there is no coupling (or no interaction between the 
ports) the S matrix of the feed network is diagonal. Therefore, 
the reflection coefficient at each port is directly the self 
reflection at that port i.e., S at 11’ and S at 22’. 

44 22 





at 


Ar\ 


When S ^ 0 and S * 0 
12 21 

11' and 22’ are quite different 
equations 4.2 and 4.3. 


, the reflection coefficient 

from S and S and are given by 
11 22 


r 

i 


s + 
1 1 


s (s r /ci -s r )) 

12 21 L 22 L ' 


...(4.2) 


r 

2 


s + 

22 


s (s r / c i-s r )) 

2i 12 a ii a 


...(4.3) 


In a two-port network, the calculation of input reflection 
coefficient at port 1, when port 2 is terminated by a load of 
reflection coefficient r , is fairly simple. Computer evaluation 

L 

is needed for a multi port network. In the present problem the 
load terminations are not simple impedances, but is an antenna 
array, which is characterized by a n-port network scatter matrix. 
The reflection coefficient at each port of the network is arrived 
at by solving a set of equations iteratively. 

Consider a 3-port network feeding a two port load as shown 

in fig 4.3. , T. refers to the reflection coefficient at port j of 
J 

network i. the primed S-parameters correspond to the second 
network and the unprimed S-parameters correspond to the first 
network. The reflection coefficient at port 3 of network 1 can be 


expressed as, 








r = S + S (a / a ) + S ( a / a ) 
a aa ai i ' a ' aa 2 a 


(4.4) 


where , 


s s r 

1 a 2i 1 


s + 

2a 


1 - 


_ 2 _ 

s r 

11 1 


1 - 2 r 


S22 + 


s s 2 r 

12 24 1 

1 - s V 

44 4 


- - • (4.5) 


and 


a 

2 r 


a 


4 

4 

s + s 

jt O 49 

2 




13 X 4 C 



a 

3 

1 - s V 

41 4 

— 

a 

3 

— 


(4.6) 


The ratios a /a , a /a are derived in appendix- I . 

2 a 1 a 

To evaluate the equation (4.4), the reflection coefficient 

of ( 2 r , Z T ) network 2 should be known. If the ports 2 and 3 of 
1 2 

network 1 are loaded by mutually Independent loads, meaning 
S-parameter of network 2 is diagonal, ( 2 r , 2 r ) are independent 

4 2 


of each other and equation (4.4) can be evaluated directly. 
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But, if the network 2 is a coupled load, (say two antenna 

elements in an array) 2 P , 2 P are also unknowns. The reflection 

a 2 

St 2 

coefficient r , are functions of S-parameters of network 2 

and reflection coefficient of network 1 , *r , 4 r , 4 r . The 

a 2 a 

expressions for 4 r , are similar to equation (4.4) to 

equation (4.6). For network-2 the reflection coefficients are. 


where , 


(b' /a' ) = S' + S (a /a' ) 

4 ' 4 ' 44 42 a 4 


• - (4.7) 


(b' /a' ) = S' + S' (a /a' ) 

2 2 22 24 4 ' 2 ' 


. . . (4.8) 


T S 

4 42 


i - s' r 

44 4 


T S 

2 24 


i - s' r 

22 2 


The reflection coefficient of network 1 4 r , 


S + S (a /a ) + S (a /a ) 

4 4 42 2 4 48 fl 4 


■ . (4.9) 


V = S + S_ ( a / a ) + S ( a / a ) 

2 22 21 i 2 22 S 2 


..(4.10) 


V = S + S (a /a ) + S (a /a ) 
a aa 34 4 3 si 2 9 ' 


. . . (4.11) 
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The incident wave ratios (a. /a.) of network 1 are functions 

»■ J 

of S-parameters of network 1 and reflection coefficients of 
network 2 (Also vice versa for (a /a.) of network 2). Hence the 

>- j 

equations to be solved are two sets of coupled simultaneous 
equations. The ratios (a /a ), (a /a ) etc., for a n-port network 

2 i 9 4. 

is lengthy, but can be calculated using a recursive relationship. 
To obtain the recursive relationship we examine the ratios (a^/a^) 
for a two port, three port and four port networks. These are 
given by the expressions derived in Appendix- I . 

Two port network : 


a 


2 


a 

1 


r 

2 


s 


21 


1 


s 

22 


r 

2 


...(4.12) 


Three port network : 



• - - (4.13) 
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Four port network : 



Where A and B are given by, 



. . . ( 4 . 16 ) 



In general , 


for an n-port network, the recursive relationship for 


30 


the ratio (a /a ) can be written as, 

m i 



r 

in 




; 1 < in 2s n 


... (4.17) 


where n represents the number of ports. 


(S. .) is given by, 
V) x 



( s J 




f s ] r s i - 

I i, x + i I I x+i , l J 

^ x+A ^ ^x-f A 

l -r f s 1 

x + i ^ x+A,x+l Jx+i — 

; x < ffin ... (4.18) 


(Si.,j) = S 

x 


x = n 


... (4.19) 



4.3. Method of Synthesis of Feed Network Characteristics: 


The configuration of the antenna system and the feed network 
are shown in fi^ 4.1. In the introduction to this chapter, we 
have discussed some preliminary points about feed network. 
Section 4.2 analyzed an n-port network, by which one can compute 
the reflection coefficient at any port, if the S-parameters 
characterizing the network and the load reflection coefficient are 
known . 


The requirement of the feed network is to provide proper 

excitations at each frequency to sustain the desired current 

distribution on the antenna. To ensure this, the driving point 

impedances and excitations have to be equal to that calculated 

earlier. In addition to this we have to ensure that the input 

reflection coefficient is minimum at the primary feed point. 

Thus, the problem Is one of minimizing the input reflection 

t. 

coefficient at (n+1) port subject to the constraints of 
providing the desired driving point Impedances and excitations to 
achieve the assumed current distribution on the antenna. 

The analysis, given so far, calculates the input reflection 
coefficient, driving point impedances and excitations, given the 
S-parameters of the antenna and the feed network. The synthesis 
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problem is an inverse problem where we have to arrive at the 
S-matrix of the feed network to give the proper driving point 
impedances and minimum reflection coefficient at the primary feed 
point. This can only be solved iteratively, using minimization 
algorithms. In this work, a pattern search method is used to 
arrive at the S-matrix of the feed network. 

In this method, an initial S-matrix is assumed for the feed 
network and its elements are successively updated based on some 
error criteria, till the error Is minimum possible. the strategy 
of pattern search[to][ n ] states that progress towards a minimum is 
accomplished b<y learning the behavior of a function in and around 
a defined base point, and based on this Information, a move is 
made. This search technique involves two type of moves viz.. 
Explore move and Pattern move. In the exploring move each 
variable of optimization is perturbed and the error is evaluated 
in the vicinity of base point. Then a vector In error space is 
generated which terminates on a point of lower error. If no lower 
error Is noted, a smaller step is used (points closer to base 
point are searched). 

The other move is pattern move. Pattern move establishes 
the new point found by the explore move and attempts to take a 
pattern step. The definition of a pattern move can be best 















Fig -d. B Flow chax-t of search tochniq'ao irn.pl intern. tod. in. DEMON. 











described by saying that it is an intelligent move based on some 
recent observations (exploratory move) of the behavior of the 
function. In a pattern move a large step in the direction of the 
vector, generated by explore move, is taken. If this is 
successful in optimizing the function, a larger step is taken and 
so on, until a failure occurs. Fig(4.4) and Fig(4.5) show the 
flow chart of the optimization algorithm. 

4. 4 Synthesis of' Feed Network Scatter Matrix: 

The formulation of [Z] matrix of the antenna array is already 


treat ed 

in 

section 3.4. 

Here 

we assume that [Z] 

matrix 

of 

the 

ant enna 

array, driving 

point 

impedances [Z ] as well 
<£& 

as 

the 

exci tat Ion 

voltages [V] 

, are 

known. As a first 

step. 

the 

[ZJ 

matrix 

Is 

converted 

to a 

corresponding [S] 

matrix 

by 

the 


transformation equations. This is a strait forward procedure 
/jfand is briefly explained here . 

[E] = [Z] [I] ...(4.20) 

where E is the voltage vector^, I is the current vector. Using 
the incident wave ’a’ and reflected wave ’b’ terminology E and I 
can be expressed as (a+b) and (a-b). Therefore, 



(a+b) 


Z (a-b) 


...(4.21) 




i£ U is the Unitary matrix then, 

S = b/a = [U+Z] _1 [Z-U] ...(4.22) 

The driving point impedance at the antenna terminals are converted 
to corresponding reflection coefficient using 


r 



...(4.23) 


Now, consider the calcu" it ion of the reflection coefficient 
at various ports of an n-port network, when it is loaded by 
another n-port network. Fig 4.6 shows the configuration used for 
calculating reflection coefficient of m th port of the n-port 
network. This involves solution of two sets of coupled 
simultaneous equations, which are given below. The subscript 1,2 
on S- parameters denote the network 1 or network 2. 

For network 1, the reflection coefficients are, 



. . .(4.24) 
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T = 
2 


21 


( 4-0 * ‘ 


s + * s 

22 23 


(±0 + ‘s ‘(-0 ♦-** s ‘(-0 

a 24 *■ a 2 n v a^ 


. .(4.25) 
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(40 
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4 r = 4 s 

n nl ^ a 


‘(- 0 *‘ s ‘(-0 ♦ ** '(-0 + *s ‘(-0 ♦... 

A J V Q, J **.«» v A <* tslA. ^ A ^ 


r»3 a 


r»4 ^ a 


+ 4 s 


For network 2 the reflection coefficients are, 


(4.27) 


2 2 
r = s 

i 




2 _ . 2 - _ 2 

+ * s U0 - z s -0 -0 
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. . .(4.28) 


T = 
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(40 
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23 


(-0 '(40 
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(40 


. . .(4.29) 
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99 
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s 

r»rt 


. . . (4.31) 


The ratios (a. / a.) can be calculated using the expressions derived 
>• i 

earlier in section 4.2, equation (4.17) to equation (4.19). 


The expressions 

(a /a.) of network 1 
' J 

_ 2 

involves 

2_ 

the 

reflection coefficients of 

network 2 ( , T » . . 

4 2 

• t 

r ), 

r» 

which 

are given by the rations 

(a /a.) of network 2. 

J 

As 

the 

ratio 

(a /a.) also involves the 

reflection coeff ici ent 

of 

network 1 


»■ J 


these equations may be solved Iteratively to obtain the reflection 
coefficient of every port of each network. The iterative process 
can be briefly described as follows. 


The reflection coefficient at port ’m' is a sum of the 
self reflection coefficients S at port m and the scatter 

mm 

parameters S . between the m th 

TYU. 


port and all the other ports 









weighted by ratios (a /a ). Hence a first order approximation can 

l W 

be made that the reflection coefficient is only due to the self 
reflection ie. , S (fig 4.7). This approximation serves as the 

lYVfTrt 

starting point for the iterative solution. The algorithm is, 

(1) Assume the load reflection coefficient of network 1 as S. 

LL 

of network 2. (This is *r , when all non dia£onal elements of 
network 2 are zero). 

(2) Solve for the reflection coefficient of all the ports of 

2 

network 1 with the terminating load reflections ( r ) as 
calculated in step (1) in first iteration and later from step (3). 

(3) Now, take the calculated reflection coefficient of network 1 
from step (2) as the terminating load reflections connected to 
network 2 and solve for the reflection coefficients of network 2. 

(4) Iterate between steps (2) and (3), till the difference 
between two successive solutions of reflection coefficient of both 
the networks are within a specified error limit. 

So far it has been discussed how to find out the 
reflection coefficient of an n-port network connected to another 
n-port network. But in our case, the antenna is an array of n 



elements characterised by n-port network Cn x n matrix!)* whereas 
the feed network has to be an Cn+1!) port network- Cn-ports 
connected to the antenna ports and Cn+ll) th port connected to the 
source!). In order to treat both the networks as a general P-port 
network* where P is the number of elements in the antenna array 
plus one* the antenna matrix can be augmented by adding an 
additional row and column with all zeros. This is equivalent to 
an Cri+l) 1 ^ port* which is uncoupled to any other ports in the 

network* and offers a matched termination. The primary feed point 
of feed network is connected to this C n +T J ’ port of antenna 

network C f i g 4.82). 

The voltage column matrix calculated in chapter -3 should be. 
next transferred in a suitable form so that it can be represented 
by the feed network £‘-par ameter s . The square of transfer 

S-par- ameter s from primary feed point* Cn+ l'J th port* to secondary 
feed points C 1 to nl> ie. * £ 2 * of the feed network 

l * <rv+-i > 

represents the power transferred from source port to antenna 

elements. From the property of Scattering parameters in a linear- 
passive non dissipative network* the power input at a port is 
equal to the sum of power transferred to all the other ports. 

Hence summing up the normalised power at secondary feed points and 
equating to one* the normalised power at primary feed point* the 
driving point voltages can be scaled down as given below: 






2 





How, the augmented Cn+1 x n+13 S-matrix of antenna is known 
and the Cn+lD th row and column of feed network is known. The rest 
of the S-parameters of feed network are given a starting value and 
the reflection coefficient at various ports [ T. to T 3 are 

t n+i 

calculated as explained earlier. 


Our object is to make these reflection coefficients looking 
into antenna ports C r to r ] equal to the corresponding 
reflection coefficient calculated earlier, which will sustain the 
necessary current distribution . To obtain a match at the primary 
feed point, the reflection coefficient at Cn+D th port should be 
zero. So, we can form the objective function as follows. 


$ 


W 

i 


2 

r 

N+i 


\ r - r 
Z 01 1 


. . . C4. 33!) 


where W^is a weighting term used to adjust the relative emphasis 
between a match to the source and a match to the driving point 
impedances. The objective function (p will be zero when the 



f eed network S-par&m^iers meet our requirements. The? ob j ecti ve 
function 4 ? is minimised using the pattern search algorithm. 

The location of antenna elements and symmetry in the 

configuration are utilised to reduce the number of variables of 

optimization. For example,, an array of seven elements as in 

Fig 4. Q* will have identical environment for elements one and 

seven, two and six, three and five. The feed network is a 

reciprocal passive network. So, S = S . From the required 

Ki ^ 

excitations at antenna terminals* S parameters are known. 

i. * rt+aL 

Wow* due to the identical array environment of different elements 
and a symmetrical current distribution about the center element* 
S = S . Considering ail these into account* the total 

*-j r,~j,rv-L 

2 

number of variables of optimisation are [ Cn+i !)/2] . 


4.5 Results and Conclusions 


This chapter started with the requirements of the feed network 
to excite an array of electrically small radiators. The problem 
is broken down to analysing an n-port network* under loaded 
conditions and then using the reverse process to synthesise the 
feed network terminal parameters. the method suggested in this 
chapter is implemented in an opt i mi s at i on routine. As an example* 
an array of five elements is designed and the results are 


presented. 


















The optimisation process started with an assumed S-matrix 
for the fee4j£wdT" network at the lowest frequency of operation. then 
the feed network S-matrix is successfully updated to meet the 
requirements based on error criteria equation £ 4 . 331 ). The 

optimisation is carried out till the error is minimum possible. 
For the next higher frequency the optimisation is repeated to 
synthesize the feed network , with the initial guess of feed 
network &-par ameter s as that of immediate lower frequency. The 
results are presented in the form of graphs £4.1 £D to £4.121). The 
£-par ameter s have a gradual variation over the bandwidth. Such a 
smooth variation is desirable from the point of view of 
realization of the feed network. the reflection coefficient at 
the primary feed point is plotted in figure £4.1^D. It is 
observed that the source remains matched over the complete 
bandwidth. The reflection coefficient at the source port is 
maximum at the lowest operating frequency. This is due to the 
hi ghl y react! ve natur e of el ectr i cal 1 y smal 1 r adi at or . £ whi eh i a 
very difficult to match as an individual radiator,, as investigated 
in section £.43. At the point when length becomes comparable to 
half wavelength, the reflection coefficient is very small 
conforming the characteristics and performance of convention 
r adi at or . 

In order to test the performance of the array and feed 
network at the frequency within the bandwidth other than one at 
which it has been designed and optimesed, the s -parameters of feed 
network, are interpolated at various points. The performance of 
the antenna are evaluated at these points and the reflection coef 
fieient at the primary feed is found to be satisfactory, 
pof ni s are also s Hown in figure £ 4 - 1 ^> • 


these 





chapter-5 Results and Conclusions 


This chapter summarises the work done and the results 
reported in this thesis and also suggestions £ or further work. 

This thesis presents a new configuration of antenna to meet 
the requirements of Direction Finding(DF) systems viz., wide band 
and constant beamwidth . The antenna is made up of electrically 
small radiators in col inear configuration, fed by a specially 
designed feed network. The feed network distributes the necessary 
excitations (as a function of frequency) in such a way that , the 
radiating currents flows over a constant electrical length 
throughout the bandwidth . 


Individual 

electrically 

smal 1 

radiator 

has 

been 

char ecteerised as a 

function of 

frequency 

using the 

well 

known 

induced EHF method . 

The results show 

that electrically 

small 


radiator is Inefficient , if operated as an individual antenna. 
Because the radiation resistance is much less compared to the 
reactance at the input terminals of the antenna. Next , mutual 
coupl ing between such electrically small radiators in colinear 
conf igurat ion is evaluated , to study the behaviour in array 
environment . In the colinear configuration, coupling is stronger 
between many adjacent elements at the lowest operating frequency 
and becomes significant only for immediate adjacent elements as 
the frequency increases . 
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An N-port network is analysed and a recursive relation is 
derived to evaluate the reflection coefficient at any port , when 
all the other (n-1) ports are loaded by a (n-1 ) port network. 
This analysis is used to synthesise the feed network terminal 
parameters . The scatter parameters of the feed network are 
arrived at using an optimisation method (pattern search) , which 
provides minimum reflections at the feed point , and excites the 
required current distribution in the array. An array of five 
elements is designed to operate over multioctave bandwidth . The 
feed network S-parameters are also synthesised for the same; and 
the results are presented in the form of graphs . The synthesis is 
carried out at discrete points over the bandwidth. By 
interpolating the synthesised S-parameters of the feed network, 
the performance of the array is evaluated over the bandwidth. The 
results show a satisfactory Input reflection coefficient over the 
entire bandwidth. 

In this work dipoles are considered as electrically small 
radiators , but this topology can be applied to any other type of 
radiating element . The design of feed network and all the 
calculations were based on [S] or [2] parameters so, the present 
work can be directly applied to any array, made up of electrically 
small radiators . Analytical formulation of self and mutual 
impedances may be difficult for complicated radiating geometry of 
individual elements In such cases, the measured scatter parameters 
can be used in the synthesis of fed network. 



The present study was restricted to equal length of elements 
of the array. One may extend this work to optimise the 
performance of the array, by choosing unequal individual elements . 
The length of the elements could be chosen to have a gradual 
reduction in radiating length. 

An important work will be to realise the feed network , 
characterized by scatter matrix. The synthesised terminal 
characteristics show a gradual variation as a function of 
frequency , and the feed network is symmetric about the center 
element. Hence it should be easy to realise the feed network. 
One may start with a conventional passive device like a T- 
junct ion or a cross (+) junction as a building block . The 
interconnection between these blocks and the charact eristics of 
the block itself can be varied to achieve the terminal parameters 
of the feed network . 
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APPENDIX i 


;CT!ON COEFFICIENT IN A GENERAL N~PORT NETWORK 


Here the ratio of incident wave amplitudes Ca^/a^3 are 
derived for a general n— port network. These ratios are evaluated for* 
3— port » 3— port and 4 —port networks > then by examining the expressions^ a 
generalised formula is developed for n-port network. 

c cause f Consider a two port network. The relationship between 

incident to reflected waves can be expressed aa » 


S a * 

A A as. 

S a + 

Si 4 


4 ^ S 


S a 
as a 


. Ca»D 
. Ca-»S3 


The reflection coefficient at port-1 C by can be expressed as. 


b 

A A 




£ Ca v'a 3 

4 2 2 A 


. C a. “"So 


Introducing the reflection coefficient of the load at port 3, 


T = Ca xb. D equation Ca-SD becomes * 
2 2 2 

b .- 'a. = S + 

£ A S:i 


C a /*’& 3 b Va. 

2 2 2 1 

Ca v*a ^ Ei^T - S 3 
2 1 2 2 2 


2A 


24 


£ Ca /a 3 

2 2 2 4 

S C a Xa "J 

22 2 4 


. „ . C a ■““40 
. Ca-B 3 
. „ . Ca-63 
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s r 

a* a 
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1 - s r 

23 2 


. C a . -“73 


Similarly aqua' 


tion C&-23 can be expressed as 

- s + 


b x a 
2 2 


22 


s Ca J 

2 1 4 2 


s r 
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where 
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c€ib& iiJ Consider* a t*hir , ee port# nstwork 
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4 A 4 
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3 4 4 
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2 2 2 
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32 2 


S a. 

4 3 3 
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- . . Ca-iOD 
. . Ca-ID 

- . . £&™i23 


3 SCI 322 BB B 

Tlie v&f 1 i on eoel fid ent# at* par t*C 13 is obbai ned f'r* om equal* i on C a— i 03 


b #»<a = 
c 4 


S 


S Ca ^a 3 + 

4 2 2 C 


S (a So. 3 


. C&dS 


4 4 i 2 2 C i 3 *3 4 

whereC are obtained from equations Ca— 113 and ta— 123 

dividing equation C a-1 23 by a 


b ^’a = S + S u .--'a 3 + S (a /a 3 . . . C a-1 S3 

a 4 ai aa a 4 as 3 4 

substituting for F = a /b in equation fa- 1B3 and bringing fa s& 3 
terms bogebher on LHS. 


l —\ ~ 
* i r . 


33 


34 
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3 2 2 4 


-163 


a F 1 
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r i 
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S + S l& s& 1 

a i-r s | 

^ Si 82 2 i J 

4 a aa 


. . . C&-173 


FrornequaiionC a-1 i 3 a^Xa.^3 eanbeobtai nedi rvthef'ol 1 owl ngmanner » 
dividing equate! on Ca— 113 by a 


b ^'a = 
a 4 


S + S € a /a 3 4- S C a .x& 3 

24 22 2 4 23 3 4 


suhsi 1 iuii m for [a ^a 1 from- equation Ca-i73 , 

si. 


. C a-1 S3 


b v”a = S +S £ a 3 +S 

2 4 24 22 2 4 23 


i - r s I 

a S3 ** 


S +S 

34 ■ 32 


Expanding t-he berm in curl, braces, 

c* c? 


3 4* S C a /a 3 + 

21 22 2 *4 


r b 5 

3 29 34 


r s 


r s s , ^ 

s 32 aa f & | 


1 - r S fa 1 

a liO V 4 •* 
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J 


C a-1 93 


Substituting for T and bringing the (a /a 3 terms to the LHS. 

2 2 4. 


f a a 1 

I — 1 

i *- J 


1 

r. 


r 3 S 3S s 23 
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The reflection coefficient at port til? can be expressed as * 
b /& ■ 5=5 S ■+• S t a /& Z? *+■ S t a /a I? ■+'23 C& /a 3 

£ ft it AS 2 £ A3 S A i4 4 A 

From equation Ca~2ED 

. . . Ca-273 

Repl aei ng b /a by t a /a 3 x’F and col 1 ecti ng the ter ms cental ni ng 

4 A *4 A *4r 

Ca /a 3 to Left Hand Side* equation (a- 273 becomes » 

<4 A 


b x’a = S + S C a ya 3 + S t a /a 3 + S t a ya 3 
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Dividing equation ta-263 by a 
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Subs t i t ut i nq f or* C a /a 3 f r om equat i on t a -223 » 
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3 1 31 32 2 4 3331 


3 4 4 i 


ta-SSU 


3* 4 


f 


1 


r s 

4 -4-* 


|S + S C a /'"a 1 + S l a --'a 3 j 

^ 42 2 i 43 3 1 j 


Ca-3CD 


Replacing b /"a by ta ^a 'J /T and collecting the terms containing 
^ 55 31 31 3 

Ca /a 5 to Left Hand Side, equation Ca-3CD becomes, 
a i 
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This completes the required formula for a four port network. 

By comparing the equations of Ca ^a D of two port network from Ca-73 
with that or 3 port network equati onC a— Si ~3 


3 f^i = 

i a i j 
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the expression for* the n-por-i network. 
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replacing S, .1 n Cn-i!) port network to S, + CS ,S f F Vti-S P I) 
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Replacing likewise for a four port network leads to the same formula 
derived equation Ca-3Q5. Hence the formula for a general n-port network 
can be written based on the observations. 


f *3_1 

i a i J 


res ] 

& 2 4 


i - r 


:s i 

22 


The ter mi inside the square bracket builds up n -times 

according to the S expression mentioned in the above paragraph. * when 

oj 


-.j 

n is the number of ports. This build up can be expressed by a mutual 
recursive f or mul aCmeani ng the same expression is substituted for each 
variable expression r* epeatedl y5 . 
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Using C&-40J Ca^/a^iof a general network can be computed. Applying the 
same logic to other- ratios Ca s /a t 5 - - - can calculate them if the 

same ratio for a network of a lower order- is known. We shall find out a 


general formula to represent the rati 


o 


Ca --'a 3 of an n-port 


network . Consider the case of Ca ^'a 3 of two, three and four- 


port 


jwtwor- ks . 

cose i .> n=2i 
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Replacing all T to T . , S to S . 

n n-KL n» n n+l ,n+l 
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Hence for n=4 we get. 
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^ 4A <42 £ 4 43 3 ± j 


Ca-4S3 


Which is same as equation ta~383 

Generalised expression can be obtai ned by combining the above? 
expression with C&-4C3. The ratio Ca /a^) of n-port network can be 

written as* 
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fa 1 
1 m I 
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l * J 


tfs p mm} m q=i 


m— 1 
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(%) 
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f a q 1 

I -CL | 

l * ± J 


. € a -483 


So for we have generalised the incident wave ratios at any 
port m with respect to incident wave at port 1 in an n-port network . 
But our interest is to calculate the reflection coefficient at any port 

l . It will be r* equi r ed to eval uaie C a Xa 1) of n -por t net wor k . Dir eci 

4 ml 

evaluation of a formula for this is complicated. As this evaluated by 
a computer* , one can simply manipulate the matrix of F and S such a way 
that the I th port is shifted to the position of por it i 3 and evaluate 
ta— 463. After the evaluation care should be taken to re-manlpulate the S 
matrix and P matrix to yield the original set. 



